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We examine the hypothesis that some supernova remnants (SNRs) may be re- 
sponsible for some unidentified 7-ray sources detected by EGRET instrument 
aboard the Compton Gamma Ray Observatory. If this is the case, 7-rays are 
CN ' produced via pion production and decay from direct inelastic collisions of accel- 

^P , erated by SNR shock wave ultrarelativistic protons with target protons of the 

interstellar medium. We develop a 3-D hydrodynamical model of SNR IC 443 
as a possible cosmic 7-ray source 2EG J0618+2234. The derived parameters of 
i-Q ' IC 443: the explosion energy Eo = 2.7 • 10^" erg, the initial hydrogen number 

*p^, density n(0) = 0.21 cm^'^, the mean radius R = 9.6 pc and the age t — 4500 yr 

O ' result in too low 7-ray flux, mainly because of the low explosion energy. There- 

^ , fore, we investigate in detail the hydrodynamics of IC 443 interaction with a 

^ ' nearby massive molecular cloud and show that the reverse shock wave consid- 

erably increases the cosmic ray density in the interaction region. Meantime, the 
Rayleigh-Taylor instability of contact discontinuity between the SNR and the 
K> , cloud provides an effective mixing of the containing cosmic ray plasma and the 

;h ' cloud material. We show that the resulting 7-ray flux is consistent with the 

5t 1 observational data. 

Key words: supernova remnants, individual: IC 443, X-rays, cosmic rays, 
J -rays 

PACS: 98.38.Mz, 98.70.Qy 98.70. Rz, 98.70.Sa 

1. Introduction 

Supernova remnants (SNRs) are believed to be the most promising accelerators 
of ultrarelativistic particles (electrons, protons and nuclei) in our Galaxy [ ||, 0. 
They are responsible for the majority of cosmic rays (CRs) with energies up to 
10^^ eV (Lorentz factor of proton 7^ ~ 10^). Classical observations of synchrotron 
radiation from SNRs, mainly in the radio band, give experimental confirmation of the 
presence of relativistic electrons in these cosmic objects [|^. Recently, the presence 
of ultrarelativistic electrons in SNRs has been confirmed by optical. X-ray and soft 
7-ray observations [ |], Q . 
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There are expected two main CR accelerators in SNR. One of them is a young 
pulsar (if there is any), i.e., a spinning neutron star with a strong magnetic field. 
The second one is a strong shock wave as the result of an interaction of supernova 
ejecta with the interstellar medium (ISM). It is the main object of our investigation. 

The theory of CR acceleration at the shock front via the first order Fermi mech- 
anism is quite well elaborated [ [^, |^, Q . Its main prediction is a considerably higher 
efficiency of the proton acceleration in comparison with an electron case, i.e., up to 
10% of kinetic energy of the incoming flow may be transformed into the energy of 
newly accelerated CRs. Experimental evidences of the proton acceleration at shock 
fronts are limited to the solar system region (a bow shock at the boundary of the 
Earth's magnetosphere, interplanetary shocks and solar flares). Therefore, it is of 
great importance to flnd additional evidences of the proton and nucleus component 
acceleration at the fronts of cosmic shock waves. 

Such a new possibility is offered by the recent observations of considerable 7-ray 
fluxes in directions on some radio bright SNRs [ §|, 0, |§] • Among the 32 unidentifled 
Galactic plane 7-ray sources detected by the Energetic Gamma-Ray Experimental 
Telescope (EGRET) aboard the Compton Gamma-Ray Observatory (CGRO) flve 
sources are coincident with SNRs without known pulsars. Four of these SNRs (W28, 
W44, 7 Cygni and IC 443) exhibit evidences of a strong interaction with interstellar 
clouds. Therefore, it is believed that the observed 7-ray emission is the result of 
an interaction of the accelerated by SNR shock wave protons with the SNR matter 
itself [ ^ or with the matter of the adjacent molecular clouds [ [7|. 

In this paper we check this hypothesis for the case of well-studied SNR IC 443 
and show that according to the atypically low initial energy of a supernova explosion 
such simple explanations do not hold here. In order to explain the 7-ray emission of 
IC 443 we elaborate a hydrodynamical model of IC 443 as a result of a supernova 
explosion in the medium with a large scale density gradient, what allows us to explain 
observational data concerning X-ray surface brightness distribution and obtain the 
parameters of remnants, especially shock wave characteristics (section 0). Using 
these parameters, we consider in detail the SNR shock interaction with a nearby 
molecular cloud (section |^) and show that the reverse shock enhances the energy and 
number density of newly accelerated CRs (section ^). Then we show that caused by 
the Rayleigh- Taylor instability, effective mixing of SNR and molecular cloud matter 
leads to the generation of an observable flux of 7-rays (section ^. Conclusions are 
given in section |^. 

2. Hydrodynamical model of SNR IC 443 

SNR IC 443 is well studied at different energy bands from radio to 7-rays [ ||, |l^, 
TT| , |12|, |T3|, [l^ (flgure [1| and 0). It is believed to be the result of a supernova explosion 



which happened about t = 3000 -i- 5000 yr ago near a dense molecular cloud at the 
distance about d Pt 1.5 kpc from us. Its unusual nonspherical morphology should 
reflect the inhomogeneity of the surrounding ISM, but until now a self-consistent 
model of IC 443 is absent. We suppose that the main reason of this is the absence 
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Figure 1. Palomar Sky Survey red plate image of IC 443. 



of a detailed numerical modelling of IC 443. So, despite the evident influence of 
inliomogeneity of the surrounding ISM, IC 443 models are mainly based on the use 
of the one-dimensional (ID) self-similar Sedov solution [ |T^ for SNR evolved in a 
uniform medium. The role of the density gradient in the surrounding ISM was taken 
into account only for the explanation of shape asymmetry within the framework of 
the approximate Kompaneetz [ |l^ approach [ 



13 



Therefore, we carry out here a detailed 2D calculation of the IC 443 evolution in 



a nonuniform medium using the proposed in [ 17, 18, 19 1 new approximate method 
of a complete hydrodynamical description of a point-like explosion in an arbitrary 
nonuniform medium. Besides geometrical sizes, we took the total flux and the surface 
brightness distribution in an X-ray band as free characteristics which give the best 
possibility for the estimation of plasma characteristics inside the SNR. 
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Figure 2. X-ray image of IC 443 obtained during the ROSAT all-sky survey [ 14| 
and 95% confidence contour (solid line) of 7-ray source 2EG J0618+2234 [ ^. 



t 



Taking into account the generally accepted distance to IC 443 d = 1.5 kpc [ 
^, [I^ , we have reproduced the following picture of the supernova explosion. At 
4500 yr ago, a SN explosion with the total energy Eo = 2.7 ■ 10^° erg took place 
in the transition zone between a warm ISM with the hydrogen number density Uo 
and an interstellar HI cloud with the number density rii and the high scale h, so the 
initial density distribution in this transition zone was n{r) = Uo + rti exp(— f//i), with 



Un 



0.16 cm 



n.: 



16 cm , h = 2.4 pc, the SN position moved on r 



13.8 



pc from the cloud centre and the initial hydrogen number density in the point of 
explosion n(0) = 0.21 cm~^. Calculations show the following best fitting parameters 
of the IC 443 model. Depending on the direction from the explosion centre, the shock 
radius R, its velocity D, the postshock temperature T and the initial hydrogen 
number density at shock front position vary in the ranges: 7.5 ^ R (pc) ^ 10.1, 
480 ^ D (km/s) ^ 880, 3.1 • 10^ ^ T (K) ^ 1.1 • 10^ 0.16 < n{R) (cm^^) ^ ^ 3^^ 
The total swept up mass of the ISM is Msnr = 28.0 masses of the Sun (M©), the 
volume of the disturbed by a shock wave region is V^ = 1.1 • 10^^ cm^, the mean SNR 
radius is i? = (3V^/47r)^/^ = 9.6 pc. In figure § the X-ray surface brightness of the 
IC 443 model is presented. The total equilibrium flux in an X-ray band is calculated 
for X-ray emissivity data from [ ^ and is obtained to be Lx{> 2.4 keV)=1.4-10^'^ 
erg/s (no absorption correction is made) which corresponds to the experimental 
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Figure 3. Surface brightness Sx distribution for the model of IC 443 in the ROSAT 
X-ray (0.1 -i- 2.4 keV) band. The angle between the symmetry axis of the SNR 
and the line of sight is 45°. The values of Sx in erg s~^ cm~^st~^ are given in 
the figure, (a) Model with no absorption correction, (b) Model which includes 
interstellar and molecular cloud absorption. 



result Lx,obs{> 2 keV) = 1.5-10^^ erg/s [ |^. In figure |Kb) we show the effect of the 
absorption of interstellar gas with the column density N^, ism = 3 • 10^^ cm"^ and 
an elongated cylindrical-like molecular cloud in front of the SNR with the radius 
Rci = 2 pc and the cloud column density Nh^ ci = (0-^5) ■ 10^^ cm~^. The data for the 
absorption cross-section from [ ^ were used. This absorption is mainly responsible 
for the unusual nonspherical X-ray shape of IC 443. 



3. Hydrodynamics of the SNR IC 443 - molecular cloud interac- 
tion 

Our calculations show that the visible X-ray shape of IC 443 is close to spherical, 
although with different radii in the south-west and the north-east, while the ROSAT 
X-ray image provides a strong evidence of two half-spheres in the IC 443 structure 



Tjl . This difference may be naturally explained by the enhanced absorption of X- 



rays in the adjacent to IC 443 elongated molecular cloud which is visible in a radio 
band [ [T( 



11, m. 22 



The same molecular cloud (or, more correctly, the part of the cloud disturbed by 
the shock wave of SNR) should be responsible for the observable 7-rays. Observations 
in CO and HCO+ emission radiolines [ ^ reveal that the total mass disturbed by 
the shock gas in the molecular cloud reaches 500-i-2000 Mq. Although this molecular 
cloud covers practically the entire remnant, the interaction region is localized in the 
southeastern part and coincides with the 7-ray source position (figure |]). It means 
that the interaction region is somewhat displaced from the SNR symmetry axis (on 
the angle -=- 30°, as viewed from the Supernova progenitor position), i.e., IC 443 
has essentially a 3D shape. 
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From the hydrodynamical point of view, the SNR - molecular cloud interaction 
is a complicated event accompanied by the arising of two shocks: a forward shock in 
the cloud matter and a reverse shock in the SNR matter. The boundary between the 
ISM and the molecular cloud becomes a contact discontinuity between the shocked 
plasmas of the SNR and the cloud [ |[ . According to the complexity of the resulting 
flow we cannot calculate it within the framework of the used above method for one 
shock flow calculations. But simple analytical estimates are possible. So, the velocity 
of the forward shock wave in the molecular cloud D^ is determined by the SNR shock 
velocity Di at the cloud position and the ratio of hydrogen number densities of the 
molecular cloud Um. and the ISM n^ [ ^: Dm = Di\/ (3ni/nm-, where (3 depends on the 
number density contrast and /5 = 6 for nm/ni — > oo. If a SNR shock wave penetrates 
into the molecular cloud at the time tp ~ 3500 yr after the SN explosion, when the 
shock radius is i?p ~ 7 pc, Di ~ 600 km/s and, according to the Sedov solution [ 
T5| , the effective thickness of the postshock gas shell is /S.Rp ~ 0.1-Rp, the typical 



time scale of the reverse shock wave existence is 

The forward shock penetrates into the molecular cloud up to the distance Lp - 
Dm ■ Atrev ~ 0-^Rp\/ Pni/um, and the shock-excited mass of the cloud is of order 



RmV ( Rp\ f Uj y/2 / Um \ V^ 

5 pc/ \7 pcj V 0.3 cm^^/ V 10^ cm~^. 



Mrn = Lp.pm-Sm-7.2.10U^] ^]--^) -^^ M,, 



where pm = ^-^nmrnp is the molecular cloud density (the typical number density 
of helium atoms riue = O.ln is taken), nip is the proton mass, Sm = t^R^ is the 
effective surface of the SNR - cloud interaction. We note that the obtained mass of 
the shocked matter of the molecular cloud is close to the observed one, if we take 
typical values n^ = 10^ -=- 10^ cm~^. 

To summarise, the proposed here 3D model of IC 443 explains in a consistent 
way the main characteristics of the remnant, including its X-ray properties and in- 
teraction with the adjacent molecular cloud. We use this model for the estimation of 
the efficiency of cosmic ray acceleration and for the explanation of 7-ray production 
in IC 443. 

4. Acceleration of cosmic rays in SNR IC 443 

CRs, i.e., high energy relativistic particles (electrons, protons, nuclei) with the 
power-law energy spectrum N oc e~°, a ~ 2.6 -r- 3.2 and the mean energy den- 
sity ujcr ~ 1 eV/cm'^ ~ 10~^^ erg/cm^ are believed to be accelerated by different 
mechanisms in our Galaxy (CRs with energies e ^ 10^^ eV) and in powerfuU ex- 
tragalactic sources (CRs with energies e ^ 10^^ eV, up to the maximal observable 
energies ^max = 3 • 10^° eV), such as active galactic nuclei, quasars, radiogalaxies 
etc. By now the most promising mechanism for the galactic CR generation is the 
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first order Fermi mechanism at shock wave fronts [ 0, ^ ^. In this mechanism fast 
particles gain energy during their diffusive motion in the shock region as a result of 
scattering on different types of magnetic field fiuctuations (magnetohydrodynami- 
cal waves, especially Alfven waves, turbulent pulsation etc.) ahead and behind the 
shock front. The main reason for effective energy gain is the convergent character of 
hydrodynamical fiows in the shock front frame. 

The main characteristics of a shock wave acceleration are the following [ [l|, 0. 
The acceleration time is 






- + - =2.1-10^7^ T7^5T--7- TTTzfr^ 77TTT-T7 U^ 



Ui — U2 \Ui U2 

where Ui {i = 1 for upstream and 2 for downstream) is the plasma velocity in the 
shock frame, ki = kc/S is the diffusion coefficient, k = rjrL^i is the mean free path, 
fL,i = s/eHi is the Larmor radius of a particle with the energy e and the charge e 
in the magnetic field Hi, c is the light velocity, D = — mi is the shock velocity in the 
observer's frame. Parameter rj is believed to be in the range 1 ^ r^ ^ 10 for SNR 
shocks. 

The theoretically predicted spectrum of accelerated particles has the power-law 
type 

N{e)de = N{eo) ■ (e/eo)""de with a = 3ui/{ui - U2) - 2, (1) 

so a = 2 for plasma with the adiabatic index 7ad = 5/3. 

Both electrons and protons are involved in the acceleration process but the 
number of accelerated electrons is considerably less than the number of protons 
Nf.{e) / Np{e) = {me/mpy~ '' . Further we shall consider only a proton component 
of SNR CRs. 

There are different estimates of the efficiency of energy transformation from the 
kinetic energy of a hydrodynamical fiow into the energy of accelerated particles [ 
|I|, ^ ^. The reasonable result is that about 10% of the kinetic energy of a fiow may 
be transformed at the shock front into the energy of CRs. The maximal energy of 
CRs is restricted by the processes of energy loss, particle escape from the region 
of acceleration, finite time of the shock existence etc. In the considered here SNR 
IC 443 case, taking into account the calculated IC 443 shock characteristics, we 
should expect the following characteristics of CR acceleration. The total energy of 
cosmic rays in IC 443 is 

„,.. . .E, ^ 10« (4) (^) e.,- 

(in order to provide the necessary amount of CRs in our Galaxy due to SNRs we 

need V = 0.03). So, for our case Eo = 2.7 • 10^° erg we obtain W^r = 8.1 • 10^^ 

For the spectral index of protons we take a = 2.1 which is close to theoretical 
predictions (|l]) and corresponds to the observable 7-ray spectrum (see below). Then 
the energy spectrum of IC 443 CRs (mainly protons) is 

„H.,.fe^(j)-°.«-:o"(5^)(^)(,,^)-" ^. ,., 
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The process of the pion creation in subrelativistic proton - proton at rest col- 
hsions is effective beginning from the proton kinetic energy e = Eq ^ 600 MeV for 
which the pion decay dominates over other mechanisms of 7-ray production [ ^. 
The same proton energy is needed to generate 7-photons with the energy e^ ^ 100 
MeV. The total number of such protons in IC 443 is 

NU> 600 MeV) = ^^ = 0.93 • 10^° (-^) (jjr^) prot. 

a — 1 Eo VO.Ol/ \10^^ eigj 

The maximal energy of CRs is restricted in the IC 443 case by the SNR age t = 4500 
yr. Then, from t = tacc we obtain 

The mean energy density of CRs inside IC 443 with volume V is 

Exactly speaking, CRs should not be uniformly distributed inside IC 443. Rather 
they will be concentrated in a thin shell of thickness AR ~ O.IR [AV ~ 0.3V) 
near the shock front, where the majority of SNR mass and about half of its thermal 



energy are concentrated [ IT5[ , since the necessary time for diffusion on the distance 



of the order of a shock radius 

'' = J^ = '°'^" (10^) (10^) ' (io5^) ^' 

considerably exceeds the IC 443 age. It means that at least half of the total energy 
of CRs Wcr is concentrated in this shell, so, the CR energy density near the shock 
front is cUgh ~ 0.5iycr/0.3V ~ l.ToJcr- As we shall see below, this value is too low for 
the generation of the necessary gamma-flux from IC 443. Therefore, we propose and 
consider in detail the promising mechanism of the CR energy density enhancement 
in IC 443, connected with the reverse shock action. 

In the reverse shock the plasma of SNR is subjected to another shock and 
its final values of the number density n/ = 7ad/(7ad — 1)'T'2 and magnetic field 



Hf = {nf /n2)H2 are 2.5 times enhanced for 7ad = 5/3 plasma [ |23|, which results 
in the increase of CR number density ncrj/ncr.2 = '^//'^2 = 2.5. Meantime the en- 
ergy of individual ultrarelativistic particles increases according to adiabatic invariant 
conservation {P]_/H = const where P± = mpV± is the normal to the magnetic field 
orientation component of particle's momentum) £f/£2 = yHf/H2. Therefore, we 
expect that behind the reverse shock front the energy density of CRs will be 



For our model we will have u^ev = 818 eV/cni'^. The maximal energy of CRs in the 
reverse shock region will be of the order of maximal energy in the main shock, i.e. 
about 10^^ eV, while acceleration conditions in the both shocks are similar. 
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5. 7-ray emission from IC 443 

The 7-ray flux from source 2EG J0618+2234 wliicli coincides witli tlie SNR 



IC 443 position is well approximated by the power-law dependence [ ^ 

-10/ ^7 ^-2-i±o.i phot 






.293 MeV/ cm2 • s ■ MeV 

^ 30 Ge^ 

kpc the 7-ray spectrum of the whole SNR is 



in the EGRET band 20 MeV ^ e^ ^ 30 GeV. For the distance to IC 443 d = 1.5 



The total 7-ray luminosity in e^ ^ 100 MeV band is L^, ohs = 2.2 • 10'^^ erg/s and 
the rate of photon creation is A^^,tot = 7.6 • lO'^^ phot/s. 

In [ P] it was shown that in the case of a pulsar absence, the most promising 
mechanism for e^ > 100 MeV 7-ray production is inelastic interaction of relativistic 
protons with protons at rest resulting in the creation of pions and their consequent 
decay into 7-rays. The 7-ray luminosity in this case is equal to the rate of energy 
transformation from relativistic protons to neutral pions [ M 



L^ = c riN / N {e)app{e)eT,-o{e)de (5) 



where n^ = lAn is the mean number density of target nuclei in the region of 
interaction, ^min ~ 600 MeV is the minimal proton kinetic energy of the effective 
pion creation (with the cross-section o-pp{e) close to the mean value a^p = 3 ■ 10^^^ 
cm^), Et^o^e) = e/6 is the mean energy transformed into the pion. 
Substituting in (0) the cosmic ray spectrum (H) we obtain: 



L^ = l^n^vW^cr = 0.2 • lO^^^n —^^ ^ . (6) 



.,_£^.,^,,.0.2.10-n(-^) ^-^ 
6 ylO*^*^ erg/ s 

So, as it was mentioned earlier in [ 0, it is necessary to have the initial number 
density in the IC 443 region n ~ 10 cm~^ and the total energy of CRs of PVcr ~ 10^° 
ergs in order to explain the observable 7-ray flux in such a simple model. But, as we 
showed in section ^ in the IC 443 case supernova exploded in a low density region 
with n ~ 0.3 cm~^ and the total energy of CRs is expected to be only of the order 
of IVcr ~ 10^^ erg. 

In [ 1^, the 7-fiux from IC 443 is explained as the result of the interaction of 
accelerated by SNR CRs with the matter of the adjacent molecular cloud. The role 
of the molecular cloud consists in the increase of the number density of target nuclei 
TT-AT in the interaction region. Namely, we can rewrite (^ in the form 

L, = ^n^iy..^^"^-^^- ^ (7) 

6 6 mp s 
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where Mgas is the total mass of the interacting with CRs gas, nip is the proton mass. 

In [ 0, this total mass is taken to be Mgas = 5 • 10^ Mq and the deduced mean 
CR energy density inside IC 443 is Ucr = 96 eV/cm^. But, as we can see from (|), 
this value of uJa is considerably higher than the one expected for IC 443. Another 
problem here is a high total mass of clouds Mgas ~ 5-10^ Mq in the region of 
interaction, while the cloud mass outside SNR is included (as we mentioned above, 
only 500 -r 2OOOM0 of gas is inside the SNR). The energy density of CRs outside 
the SNR should be considerably lower than that inside it [ |I], §], and, therefore, the 
pion production should be depressed. 

We suppose that the observable 7-ray flux from IC 443 is caused by enhancement 
of the CR energy density in the reverse shock from the SNR - molecular cloud 
interaction. Indeed, using (^ we obtain from (0) the following 7-ray luminosity in 
the case of enhancement of the CR energy density: 



L^ = 6.8 • 10-2^^ P M,^W,,R^ 
nip V^2/ ' 



3/2 

3 
p 



2 9 • 10^^ ' W^r \ f M^ \ f Rp\ ^ erg 



1048 ergy/ \^103 MqJ V7pc 

Therefore, we obtain the theoretical 7-ray luminosity Lj = 2.3 • 10^^ erg/s for IC 443 
despite the low explosion energy Eo = 2.7 • 10^° erg and, accordingly, the low total 
energy of cosmic rays Wcr = VEo = 8.1 • 10^^ erg. 

Additional support for this model follows from the existence of an effective mix- 
ing mechanism of a hot containing CR plasma and dense shocked plasma of the 
molecular cloud - the well known Ray leigh- Taylor instability of a contact disconti- 
nuity between the both media which should lead to a decay of contact discontinuity 
and fragmentation of the shocked cloud material onto a set of cloudlets [ |T^ . 

Another possibility of the effective CR - cloud material interaction takes place 
in the case of a strongly nonuniform cloud, in which dense molecular cores of mass 
~ 10 Mq and number density ~ 10^ cm~^ are embedded in a relatively less dense 
intercloud medium with the number density ~ 10 ^ 10^ cm~^. As a result of dis- 
placement of a contact discontinuity in the molecular cloud direction, such separate 
cores penetrate into the shocked by a reverse shock plasma region and should be the 
sites of an intense 7-ray production. Radio observations of the shocked molecular 
gas in the IC 443 region reveal the existence of few such cloudlets [ 52 . 



6. Conclusions 

If the association of some unidentified EGRET sources with SNRs is confirmed, 
it will be a new argument for the CR acceleration in SNRs. It is of great importance 
that contrary to a radio band, 7-rays give unique information about the acceler- 
ation of a proton component of CRs in SNRs. In this work we have proposed a 
self-consistent model of IC 443 which explains a total set of its main parameters, 
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including shape, size, X-ray and 7-ray fluxes. We have shown that, despite the low 
supernova explosion energy, account for the SNR - molecular cloud interaction, es- 
pecially for the role of the reverse shock wave, is the main reason of an intense 7-ray 
flux from IC 443. Further observational and theoretical efforts are needed for check- 
ing the theory of CRs acceleration by SNRs and IC 443 is a prominent candidate 
for such investigations. 
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